Study of the process e + e — ► LJ7T with the 

KLOE detector 



The KLOE Collaboration 

F. Ambrosino d , A. Antonelli a , M. Antonelli a , F. Archilli a , 
C. Bacci h , P. Beltrame b , G. Bencivenni a , S. Bertolucci a , 
C. Bini g , C. Bloise a , S. Bocchetta h , V. Bocci g , F. Bossi a , 
P. Branchini h , R. Caloi g , P. Campana a , G. Capon a , 
T. Capussela d , F. Ceradini 11 , S. Chi a , G. Chiefari d , 
P. Ciambrone a , E. De Lucia a , A. De Santis g , P. De Simone a , 
G. De Zorzi g , A. Denig b , A. Di Domenico g , C. Di Donato d , 
S. Di Falco 6 , B. Di Micco h , A. Doria d , M. Dreucci a , G. Felici £ 

A. Ferrari a , M. L. Ferrer a , G. Finocchiaro a , S. Fiore g , 
C. Forti a , P. Franzini 8 , C. Gatti a , P. Gauzzi g , S. Giovannella £ 
E. Gorini c , E. Graziani h , M. Incagli 6 , W. Kluge b , V. Kulikov k 

F. Lacava g , G. Lanfranchi a , J. Lee-Franzini a,1 5 D. Leone b , 
M. Martini a , P. Massarotti d , W. Mei a , S. Meola d , S. Miscetti a 
M. Moulson a , S. Miiller a , F. Murtas a , M. Napolitano d , 
F. Nguyen h , M. Palutan a , E. Pasqualucci g , A. Passed h , 
V. Patera a ' f , F. Perfetto d , M. Primavera c , P. Santangelo a , 
G. Saracino d , B. Sciascia a , A. Sciubba a ' f , F. Scuri e , I. Sfiligoi a 
T. Spadaro a , M. Testa 8 ', L. Tortora 11 , P. Valente g , 
B. Valeriani b , G. Venanzoni a , R.Versaci a , G. Xu a 'j 

a Laboratori Nazionali di Frascati delVINFN, Frascati, Italy. 

h Institut fur Experimentelle Kernphysik, Universitdt Karlsruhe, Germany. 

c Dipartimento di Fisica dell'Universita e Sezione IN FN, Lecce, Italy. 

d Dipartimento di Scienze Fisiche dell'Universita "Federico II" e Sezione INFN, 

Napoli, Italy 

c Dipartimento di Fisica dell'Universita e Sezione INFN, Pisa, Italy. 

* Dipartimento di Energetica dell'Universita "La Sapienza" , Roma, Italy. 

g Dipartimento di Fisica dell'Universita "La Sapienza" e Sezione INFN, Roma, 

Italy. 

h Dipartimento di Fisica dell'Universita "Roma Tre" e Sezione INFN, Roma, Italy. 



1 Physics Department, State University of New York at Stony Brook, USA. 

3 Permanent address: Institute of High Energy Physics of Academica Sinica, 

Beijing, China. 

^Permanent address: Institute for Theoretical and Experimental Physics, Moscow, 

Russia. 



Abstract 

Using ~ 600 ptr 1 collected with the KLOE detector at DA$NE, we have studied 
the production cross section of 7r + 7r~7r°7r° and 7r°7r°7 final states in e + e~ collisions 
at center of mass energies between 1000 and 1030 MeV. By fitting the observed in- 
terference pattern around for both final states, we extract a measurement (pre- 
liminary) for the ratio T(u -► 7r°j)/T(uj -► vr+vr-vr ) = 0.0934 ± 0.0022. Since these 
two final states represent the 98% of the to decay channels, we use unitarity to derive 
BR(w -► tt+tt-tt ) = (89.94±0.23)% and BR(w -► tt°7) = (8.40±0.19)%. Moreover, 
the parameters describing the e + e~ — > 7r + 7r~7r°7r° reaction around M<p are used to 
extract the branching fraction for the OZI and G-parity violating 0^u;7r o decay: 
BR(0 -► lj7t°) = (5.63 ± 0.70) x 10" 5 . 

Key words: e + e~ collisions, rare (j) decays, VMD, OZI violation, Isospin violation 



1 Introduction 

In the energy region of few tens of MeV around M^, the e + e" — > 7r + 7r~7r°7r° 
production cross section is largely dominated by the non-resonant processes 
e + e~ — > p/p' — > ujn°. However, in a region closer to M^, a contribution from 
the OZI and G-parity violating decay (fi^um is expected. This strongly sup- 
pressed decay can be observed only through the interference pattern with the 
previous reaction, which shows up as a dip in the production cross section as 
a function of the center of mass energy (\/s). A fit to the cross section energy 
dependence allows us to extract the (fi^un branching fraction (BR). 

There is a much more complicated interference scenario for the e + e~ — > 7r°7r°7 
channel. Here we expect contributions also from <fi — > pi\ and — > S7 inter- 
mediate states, where S is a scalar meson. In another paper [1] we have shown 
that at yfs ~ the interference between — > Sj and e + e~ — > ujti events, 
evaluated by fitting the M^-M^ Dalitz plot, is small. Assuming this effect to 
be negligible to first order, a fit to the cross section interference pattern for this 
final state will nevertheless give information about the e + e~ — > p/p' — > ujtt 
process and the resonant decays <fi — > loh® and <p — > p7T°. Comparing with the 
fit to the 7r + 7r~7r°7r° channel, the ratio T{u) — > 7r°7)/T(u; — > 7r + 7r~7r°) can be 
extracted. 
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2 The KLOE detector 



The KLOE [2] experiment operates at DA$NE [3], the Frascati 0-factory. 
DA$NE is an e + e~ collider running at a center of mass energy of ~ 1020 MeV, 
the mass of the 0-meson. Equal-energy positron and electron beams collide at 
an angle of 7T-25 mrad, producing 0- mesons nearly at rest. 

The KLOE detector consists of a large cylindrical drift chamber, DC, sur- 
rounded by a lead-scintillating fiber electromagnetic calorimeter, EMC. A su- 
perconducting coil around the EMC provides a 0.52 T field. The drift cham- 
ber [4], 4 m in diameter and 3.3 m long, has 12,582 all-stereo tungsten sense 
wires and 37,746 aluminium field wires. The chamber shell is made of car- 
bon fiber-epoxy composite and the gas used is a 90% helium, 10% isobu- 
tane mixture. These features maximize transparency to photons and reduce 
K L — > K s regeneration and multiple scattering. The position resolutions are 
o"j; y ~150 fim and <r 2 ~ 2 mm. The momentum resolution is cr(p±)/p± ~ 0.4%. 
Vertices are reconstructed with a spatial resolution of ~3 mm. The calorime- 
ter [5] is divided into a barrel and two endcaps, for a total of 88 modules, 
and covers 98% of the solid angle. The modules are read out at both ends 
by photo-multipliers, both in amplitude and time. The readout granularity is 
~ (4.4x4.4) cm 2 , for a total of 2440 cells . The energy deposits are obtained 
from the signals amplitude while the arrival times of particles and the posi- 
tions in three dimensions are obtained from the time differences. Cells close in 
time and space are grouped into a calorimeter cluster. The cluster energy E 
is the sum of the cell energies. The cluster time T and position R are energy 



weighed averages. Energy and time resolutions are &e/E = 5.7%/ J E (GeV) 



and a t = 57 ps/^E (GeV) © 50 ps, respectively. The KLOE trigger [6] uses 
both calorimeter and chamber information. In this analysis all the events are 
selected by the trigger calorimeter, requiring two energy deposits with E > 50 
MeV for the barrel and E > 150 MeV for the endcaps. A cosmic veto reject 
events where at least two outermost EMC layers are fired. 



3 The ^fs dependence of e + e — > 7T + 7T 7r 7r°/7r 7r°7 cross sections 

As mentioned before, in the energy region below 1.4 GeV the ii + it~it it /it it ^ 
production cross sections are dominated by the non-resonant process 
e + e~ — > p/p' — > um°. At \/s ~ M^, the decay — > un also contributes 
and interferes with the other processes. In the neutral channel there are also 
contributions from — > S7 and — > pn°. The dependence of the cross section 
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on the center of mass energy can be parametrized in the form [8]: 



(t(Vs) = <To(y/s) 



1 - Z- 



where o"o(\/i) is the bare cross section for the non- resonant process, Z is the 
complex interference parameter (i.e. the ratio between the decay amplitude 
and the non resonant processes), while and D<p are the mass, the width 

and the inverse propagator of the meson respectively. The non-resonant cross 
section in this energy range increases linearly with y/s. A model independent 
parametrization 2will be used in this paper. 

<To(Vs) = a + a'(^- M^) (2) 



4 Data analysis 



All the available statistics collected at the peak in 2001-2002 data-taking 
periods, corresponding to 450 pb _1 , has been analyzed. Moreover four scan 
points (1010 MeV, 1018 MeV, 1023 MeV and 1030 MeV) of ~ 10 pb" 1 each 
and the off-peak data (1000 MeV) acquired in 2005-2006 have been included 
in this analysis. All runs are grouped in center of mass energy bins of 100 keV. 

4-1 e + e~ — > um° — > 7r + 7r~7r°7r° 



In the 7r + 7r~7r°7r° analysis, data are filtered by selecting events with the ex- 
pected final state signature: two tracks connected to a vertex inside a small 
cylindrical fiducial volume around the Interaction Point (IP) and four neu- 
tral clusters in the prompt Time Window (TW), defined as |T 7 — R y /c\ < 
MIN(3.5 ■ Gti 2 ns). To minimize contamination from machine background, 
we require also a minimal energy for the clusters (10 MeV) and a minimal 
angle with respect to the beam line (~23°). A kinematic fit requiring total 

Table 1 

Background channels for W7r ^7r + 7r~7r 7r . Signal over background ratios after ac- 
ceptance selection and xicfit cut are reported for events collected at ^fs ~ M^. 



Channel 


S/B (acc) 


S/B (x! fit cut) 


K S K L 


1 


10 


K+K- 


10 


60 


pTT 


30 


200 




20 


800 
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four-momentum conservation and time of flight (TOF) for photons improves 
the energy resolution. The resulting \ 2 (xk&t) * s use d to select a signal enriched 
(Xxfit < 50), S cvts , and a background dominated (xkfit > 50), B cvts , samples. 
The signal analysis efficiency in the S evts sample has been evaluted by Mon- 
tecarlo (MC). The resulting value e ~ 40% is dominated by the acceptance 
requirements and has a small dependence as a function of y/s. 

The background channels are listed in Tab. 1. The main contribution comes 
from <fi — > K S K L — > 7r + 7r~7r 7r° and — > K + K~ with — > 7r ± 7r°, which 
have the same final state. The first one has also a comparable production cross 
section with respect to the signal at the <fi peak. The other two background 
components ((f) — > rj'-f with rj — > 7r + 7r~7r°, and — > 7r + 7r _ 7r°) mimic the final 
state signature because of additional clusters due to accidental coincidence of 
machine background events and/or shower fragments (splitting). In the signal 
enriched region, the expected contamination at y/s ~ is ~ 12%. 



1 

& 60000 



40000 



20000 



> 

o 



M 30000 - 



a 

> 20000 



10000 - 




M rec (MeV) 

Fig. 1. Data-MC comparison for 7r + 7r~7r°7r° signal enriched distribution using events 
taken at 1019.75 MeV : (a) XKfit (Ndf=8); (b) cosine of the angle between recon- 
structed 7r°'s; (c) 7T° recoil mass. Black dots are data, while hatched and white 
histograms are MC signal and background shapes respectively, weighted by our fit 
results. 



The signal counting on data is performed for each y/s bin by fitting the re- 
constructed 7T° recoil mass (M rec ) distribution for both S evts and B evts samples 
with MC signal and background shapes. The fit procedure is based on a likeli- 
hood function which takes into account both data and MC statistics. In Fig. 1, 
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data-MC comparison of few relevant distributions for the most populated en- 
ergy bin is shown. 

The results are summarized in Tab. 2 where the signal counts, the x 2 of the fit 
and the visible cross section are reported for all \fs bins. A preliminary esti- 
mate of the systematic error, dominated by tracking and vertexing efficiency, 
has been added to the a% error. 

Table 2 

Signal counting, x 2 of the fit and visible cross section for e + e~ — > tt + tt^tt°tt events. 
The errors on er^ contains a relative systematic error contribution of 1.8%. 
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1.05 


5.81 ±0.11 
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1.10 


5.74 ±0.11 
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336385 ± 2271 


1.04 


5.86 ±0.11 


1019.85 


264184 ± 2061 


1.06 


6.05 ±0.12 


1019.95 
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1.00 


6.16 ±0.15 


1020.05 


18358 ± 433 


1.01 


6.10 ±0.18 


1020.15 


7293 ± 356 


0.99 


6.06 ± 0.32 


1020.41 


9067 ± 222 


0.92 


6.08 ±0.19 


1022.09 


19307 ± 242 


1.06 


7.08 ±0.16 


1022.98 


29995 ± 265 


0.93 


7.52 ±0.16 


1029.97 


35125 ± 489 


1.01 


8.03 ±0.22 



4.2 e+e- -> U7i° -> 7r°7r°7 

The acceptance selection for 7r°7r°7 events requires five neutral clusters with 
E 1 > 7 MeV and a polar angle |cos# 7 | < 0.92 in the prompt Time Win- 
dow. After applying a first kinematic fit (Fitl) imposing total 4- momentum 
conservation, photons are paired to 7r°'s, by minimizing a x 2 built using the 
invariant mass of the two 77 pairs. A second kinematic fit (Fit2) imposes also 
constraints on the 7r° masses. 

The background with final state different from 7r°7r°7 is rejected by requiring 
x| it2 /Ndf < 5 and AM 77 = |M 77 — M n \ < 5(T 77 , where M 77 and <r 77 are 
evaluated using the photon momenta from Fitl. After these cuts the remaining 
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sample is dominated by e + e~ — > luti — > 7r°7r°7 and — > S — > 7r 7r°7 events. 
Signal is then selected neglecting the interference between the two processes 
and cutting on the intermediate state mass. Since M 7r7 the closest mass to 
M u of the two 7r°7 combinations, only events satisfying \M nj — M u \ < Sctm^ 
are retained. In Tab. 3, the background channels are listed together with the 
S/B ratio before and after the application of the whole analysis chain. The 
residual 10% background contamination comes predominantly from <fi — > 777 — > 
7r°7r 7r°7 events where two photons are lost or merged. 

In Fig. 2 data-MC comparison for events in the most populated y/s bin is 
shown. The ip variable is the minimum angle between the photon and the 7r°'s 
in the di-pion rest frame. A good agreement is observed both after acceptance 
selection and after applying analysis cuts. 
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Fig. 2. Data-MC comparison for events taken at 1019.75 MeV. Top: normal- 

ized x 2 °f the second kinematic fit after acceptance cuts. Bottom: 7r°7 invariant mass 
(left), and cos rp distribution after cutting on M n ^ (right). In the upper panel all the 
background is grouped together, while in the lower ones the (j) — > 5*7 contribution 
is shown alone. 
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The overall analysis efficiency for the identification of the signal is evaluated 
by applying the whole analysis chain to signal MC events: e nnl ~ 40%, almost 
flat in \/s. The value obtained for each bin, together with the corresponding 
integrated luminosity, has been applied to the signal counting to obtain the 
visible cross section. Results are summarized in Tab. 4; errors include statistics 
and background subtraction only. 



Table 3 

Background channels for e + e" The signal over background ratio 
before and after the application of the analysis cuts is reported for events collected 
at y/s ~ M(j>. 

Background S/B (no cuts) S/B (all cuts) 



4> -► Sj ->■ 7r°^°7 1.5 35 

cj) — ► r/7r°7 — ► 777r°7 5.4 120 

-► 7/7 -► 7r 7r°7r°7 0.04 15 

-► 7/7 -► 777 0.04 380 

-► vr°7 0.13 2840 



Table 4 

Signal counting and visible cross section for e + e~ — ► 7r°7r°7 events. 
(MeV) N^t ± <5jv <T±^ 



1000. 


10 


5523 ± 75 


0. 


540 


± 


0. 


.007 


1009. 


.90 


2445 ± 50 


0. 


.607 


± 


0. 


.012 


1017. 


.15 


831 ± 30 


0. 


.593 


± 


0. 


.020 


1017. 


,25 


680 ± 28 





.578 


± 


0. 


.022 


1018. 


20 


2088 ± 50 


0. 


.554 


± 


0. 


013 


1019. 


.35 


273 ± 18 





.547 


± 


0. 


036 


1019 


.45 


4911 ± 79 


0. 


.514 


± 


0. 


.008 


1019. 


.55 


7693 ± 100 





.510 


± 


0. 


.006 


1019 


.65 


14788 ± 141 


0. 


.532 


± 


0. 


.005 


1019 


.75 


27556 ± 199 


0. 


.530 


± 


0. 


.004 


1019 


.85 


20927 ± 170 


0. 


.529 


± 


0. 


.004 


1019 


.95 


2869 ± 60 


0. 


.528 


± 


0. 


Oil 


1020 


.05 


1475 ± 43 


0. 


.536 


± 


0. 


.015 


1020 


.15 


577 ± 26 





.536 


± 


0. 


.024 


1020 


.45 


542 ± 26 


0. 


.524 


± 


0. 


.024 


1022 


.25 


996 ± 33 


0. 


.639 


± 


0. 


.021 


1022 


.35 


415 ±21 


0. 


.661 


± 


0. 


.034 


1022 


.95 


2574 ± 53 





.646 


± 


0. 


013 


1029. 


.95 


3233 ± 57 





.751 


± 


0. 
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5 Fit results and uj branching ratios extraction 



The measured values of visible cross section, shown in Tab. 2 and 4, are fitted 
with the parametrization (1), convoluted with a radiator function [7]. The free 
parameters are: a l Q , §t(Zi), ^s(Zi) and a[, where % is the 47r or nn^ final state. 
In Fig. 3 data points with the superimposed fit function are shown for both 
channels. The preliminary values for the extracted parameters are reported in 
Tabs. 5. The resulting x /Ndf are 12.8/15 (P(x 2 ) = 62%) for the fully neutral 
channel and 13.4/13 (P(x 2 ) = 42%) for the other one. 
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Fig. 3. Cross section fit results for the e + e — > ir + ir 7r°7r° (top) and e + e 
(bottom) channels. Black dots are data, solid line is the resulting fit function. 
Table 5 

Fit results for the e + e~ — > tt + tt~tt°tt cross section (left) and for e + e~ 
cross section (right). 



7r°7r° 7 



^o 7r o 7 



Parameter (e + e — > ir + ir 7r°7r°) 

*t (nb) 
»(Z 47r ) 

9(Z 47r ) 

<r 47r (nb/MeV) 



Parameter (e + e — ► 7r°7r°7) 



8.12 ± 0.14 
0.097 ±0.012 
-0.133 ±0.009 
0.072 ± 0.008 



7T7T7 



(nb) 



^(■^-777) 

a' nni (nb/MeV) 0.0079 ± 0.0006 



0.776 ± 0.012 
0.013 ±0.013 
-0.155 ±0.007 



From the two previous measurements we obtain: 



a (uj -> tt°7) 
cr (u; — > 7r + 7r~7r°) 



0.0956 ±0.0022 



(3) 
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Taking into account the phase space difference between the two decays [8] , the 
ratio of the partial widths can be extracted: 

— [; = 0,0934 ±0.0021 (4) 



Since these two final states the 98% of the to decay channels, we use the 
T(uj — > 7T°7) /T{u — > 7T + 7r _ 7r ) ratio and the sum of averages of the existing BR 
measurements on rarest decays [9] to extract the main uj branching fractions, 
imposing the unitarity relation: 



BR{u -> 7T+7r-7r ) = (89.94 ±0.23)% (5) 
BR(oj^7t 7) = (8.40 ±0.19)% (6) 



with a correlation of 82%. Comparison between our evaluation and the values 
in PDG [9] is shown in Fig. 4. 

.-90.5 i 1 




8 8 - 5 BR(CQ 9 -WV) (%) 

Fig. 4. Branching fraction for the two main u decay channels. The black square 
is the KLOE fit result, while the black dot is the PDG constrained fit result. The 
shaded regions are the 68% C.L. 
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6 BR(0 — > c<j7r ) evaluation 

The measured and Z 47r paramters of the 7r + 7r~7r°7r° final state are related 
to the BR(0^cj7t°) through the relation: 

BR{<P - cvr ) = ffoK)|Z4 - |2 ; (7) 

where cro(m^,) is the total cross section of the e + e _ — > c<j7r process and is 
the peak value of the production cross section for the <f> resonance. 

Using the parameters obtained from the 7r + 7r~7r°7r° analysis, the T ee measure- 
ment from KLOE [10] for the evaluation of o^, and our value for BR{uo — > 
7r + 7r~7r°) we extract: 

BR((f) -> UTv°) = (5.63 ± 0.70) x 10~ 5 (8) 
in agreement with the previous measurement from the SND experiment [8]. 
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